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Abstract 



1 Introduction 



As suggested by London ||T| and substantiated by Feynman there is an intimate 
relation between the lambda transition in ^He and the Bose-Einstein condensation 
of the ideal Bose gas (IBG). The IBG explains some basic properties of liquid 
^He like for example the irrotational superfluid flow. We propose a modification 
of the IBG that leads to a realistic expression for the specific heat. 
The single particle functions of the IBG are of the form 

ipk oc exp(ik ■ r) = exp[i {kix + + k^z)]. (1) 

Without changing the kinetic energy we may alternatively consider the single par- 
ticle functions 

OC exp [i{k2y + ksz)] sm{qx + 0), (2) 

where q = ki. For an arbitrary direction (chosen here as the x direction) the 
functions ipl^°' have a specific phase 0; they are phase ordered (p.o.). This phase 
ordering implies the following correlation: For ki = k[ = q, the probability densi- 
ties Iv^k" P ^^"^ lv^k'° P have common minima and maxima. That means that two 
particles with the same momentum in one extinguished direction have an addi- 
tional positive spatial correlation. No such correlations are present for the single 
particle functions (|I|). 

The phase ordering leads to a correlation energy of the form 

-Ecorr OC ^ with Ug = ^ U]^. (3) 

Here is the number of atoms with momentum k = (5,^25^3) and Ug is the 
number of atoms with common momentum g in a; direction. In the group of Ug 
atoms each one has an extra correlation with the other atoms; this results in a Ug 
contribution to the energy. This correlation energy has been investigated in detail 
in Ref. §. 

We will use the IGB many-body wave functions together with Jastrow factors 
as proposed by Chester [Q. The effective interaction on the model level is then 
attractive (the hard core is cut out by the Jastrow factors). This means that the 
energy -Ecorr is negative and that the correlations are favoured by the condition of 
minimal free energy. 

For a finite correlation energy one must assume that the phase ordering is local. 
The phase correlations of the single particle functions @ could be enforced by 
physical boundary conditions for the considered macroscopic volume V (instead 
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of periodic boundary conditions for the single particle functions (|1|)). In this 
case the correlations are a surface effect implying ii^corr oc V^^/^, and Ecow/V 
for y — > oo. Assuming local phase ordering introduces a parameter Vq that 
defines the finite range over which the single particle functions are correlated. The 
directions (x direction in Eq. (0)) of phase ordering are in general different in 
various subvolumes. The statistical average over these subvolumes ensures the 
isotropy of the bulk liquid. Evaluating the energy with the IBG expectation 
values (nk) yields a logarithmic singularity 0, i.e., (-Ecorr) — )f:ln|)f:|, where t 
is the relative temperature. Adjusting the amplitude of this singularity to the 
experimental value fixes the model parameter Vq (leading to Vq/v ^ 10^, where v 
is the volume per atom). 

It is the aim of this paper to calculate the entropy change due to the assumed 
local phase ordering. This calculation is independent of the previous determination 
of the correlation energy (|^). In order to calculate the phase ordering entropy we 
establish a relation between the functions and coherent states; this relation 
is analogous to Anderson's connection |^ between the condensate wave function 
and locally coherent states. For coherent states the phase variances are given 
by A0k = l/(2A/nk). Restricting the phases by A0k yields the phase ordering 
entropy S'p,o(nk). Using the occupation numbers (nk) of the IBG form we find 
logarithmic singularity for the thermodynamic entropy, (5'p.o) oc —t In \t\. For the 
specific heat per particle, c = —A In we obtain the result A = 3k^/2, where 
is Boltzmann's constant. 

As a theory of the lambda transition, the presented approach is a modest 
attempt only. Instead of determining exact many-body states we make a guess of 
the relevant correlations. Instead of evaluating the partition sum we use the IBG 
expectation values (arguing that the correlation effects are of minor influence on 
the occupation of the individual levels). The determination of the exact many- 
body states and the exact evaluation of the partition sum is, of course, not feasible. 
One might, however, try to find a suitable variational ansatz for the many-body 
states for which the free energy can be calculated and minimized (like it may 
be done for BCS states). The present evaluation of the entropy together with 
the previous evaluation of the correlation energy (^ are a first step towards such 
an approach. As it stands, we present a phenomenological many-body model that 
yields a promising result, namely a logarithmic singularity with a sensible strength. 

As a many-body approach for the lambda transition the presented model lies 
outside the main stream of realistic models for liquid helium. There are sev- 
eral well-known and thoroughly studied microscopic many-body approaches (for 
an overview see Ref. 0), hke diagrammatic (Green's function) approaches [|], 
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variational methods (in particular the correlated basis function method and its 



generalizations ||^, 10|), and Monte Carlo calculations. As far as the approaches 
are based on elementary excitations they are restricted to low temperatures. For 
temperatures around the lambda point, there are successful applications of path- 
integral Monte Carlo simulations [jTT], and of the variational approach [|l^]. 



These approaches do, however, not yield analytic and realistic expressions for the 
asymptotic behavior. What is nowadays widely considered as the theory of the 
lambda transition and its asymptotic properties is the renormalization group the- 
ory [Q. In contrast to the model presented here, this theory is not a many-body 
description. 

Several authors [l^, |18[ have considered coherent states in the descrip- 

tion of liquid helium, too. Their investigations concern mostly the connection 
between the condensate and coherent states as in Anderson's paper p. The local 
phase ordering that we consider is not specifically related to any of these other 
investigations. 



2 Phase ordering entropy 
2.1 Many-body states 

Chester |^ studied the density matrix for many-body wave functions of the form 
\1/ = 5F\l/iBG- Here S denotes the symmetrization operator, F = Ylij fi.'fij) is the 
Jastrow factor, \l/iBG(rj, ""-k) describes the IBG state, and defines the positions 
of the helium atoms. Using hard sphere factors for /(r), Chester obtained the 



IBG values (r;,k)iBG for the average occupation numbers (rik). McMillan used 
the ansatz \1/ = iSF for the ground state; for realistic atom- atom interactions he 
determined the function /(r) by minimizing the energy. This yields a realistic pair 
correlation function and a realistic binding energy. 

We follow basically Chester's approach replacing, however, the single particle 
functions (|1|) by the phase ordered ones (0). The many-body wave function is then 
of the form 

N 

'^ = SFX{vi:-{v,)=SFX[[vrr. (4) 

3=1 k 

The factor [v^k]'^'' in the last expression stands for the product V9k(ri^+i) ■ (y3k(ri^+2) ■ 
. . . ■ V3k(rj/+nfc)- The total number of atoms is denoted by A^. 

As explained in the Introduction, a finite correlation energy (^) requires finite 
extensions of the phase ordered single particle functions. For this purpose the 
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volume V is thought to be divided into V/Vq subvolumes of size Vq, and the single 
particle functions are restricted to such subvolumes. The phase ordering can then 
be obtained by requiring physical boundary conditions for the subvolumes, or by 
considering coherent states in these subvolumes (Sec. pl3|) . The product over the 
subvolumes might be explicitly displayed in the last expression of Eq. (^) . We will 
adjust the prefactor of the momentum integrals such that the result corresponds 
to the total volume V. 

The Jastrow factors are of no influence on the statistical counting. Without 
phase ordering (i.e., using the single particle functions (|l]) in Eq. (||)) one obtains, 
therefore, the well-known IBG expression for the entropy: 

•S'iGB = ^ [(1 + nk) ln(l + nk) - Uk In(nk)] . (5) 

k 

The entropy change due to the phase ordering may be divided into two parts: 

(a) All Ug atoms in states with ki = q must adopt the same (mean) phase. This 
leads to a replacement of the entropy (|]) by a modified expression S[-qq. 

(b) A well-defined phase requires a small phase variance. The entropy change 
due to this restriction is denoted by Sp,o.- This decisive contribution is called 
phase ordering entropy. 

2.2 Mean phase 

Admitting the values = and = 7r/2 in Eq. (0) means that we use the functions 
sin(gx) and cos(ga;). This is equivalent to the use of exp(ifcia;) and exp(— ifcix), 
i.e., to the case (0) of no phase ordering. Using the single particle functions (|^) 
with one definite phase (lets say = 0) effectively implies that only every second 
single particle state is occupied. This can be accounted for by the substitutions 
~^ (1/2) J2k '^k ~^ in the expression (^, i.e., 

Smc = Y 5Z + 2^k) ln(l + 2n^) - 2nk ln(2nk)] . (6) 

k 

The step from the entropy (H) to the entropy @ takes into account the choice 
= for all atoms. Deviating from this, we could admit different 0^ values in 
Eq. (^ for different g's without destroying the considered correlations. This is, 
however, a negligible effect because the increase in entropy would be AS/N = 
Eg---/N ^ 0(Ar-2/3) ~ 0. 

The expectation value of S{^q has similar temperature dependence as that of 
^iBG (Sec. 13). 
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2.3 Phase variances 



Anderson has established a connection between the condensate wave function 
and locally coherent states. The coherent states are constructed in suitably chosen 
subvolumes AV^. The volumes must be large enough in order to define a mean 
phase (f) with a small variance A(f) ^ 1; this requires that the number of atoms in 
AV is large compared to 1. On the other hand, the volumes AV must be small 
enough in order to ensure ~ const, in spite of macroscopic perturbations or 
flows. 

The connection between the condensate wave function and locally coherent 
states may be carried over to noncondensed phase ordered states (@) provided that 
nk ^ 1. This means that we may construct coherent states that correspond to 
the single particle functions (|]). As in Anderson's work these coherent states are 
restricted to flnite volumes. For the following calculation it is not necessary to 
quantify the size of these subvolumes. 

Appendix presents an explicit construction of coherent states that correspond 
to phase ordered single particle functions @. Here we restrict ourselves to the 
presented argument based on the analogy to known work 0. 

If the atoms with momentum k form a coherent state then their phase variance 



20|| reads 



A0k = -j= . (7) 

The bar denotes the quantum mechanical average. 

In Eq. (|T]) we may admit arbitrary phases (pj for each atom, v^k(i"j) oc exp(ik- 
Tj + (pj). These phases 4>j are of no influence on physical quantities (like the single 
particle kinetic energy or an interaction matrix element). In particular, they are 
of no influence on the statistical counting. Such phases are usually ignored. 

Contrary to this, arbitrary phases (pj in Eq. (|^), {p^°' oc sm{qxj + would 
destroy the considered correlations. The phase ordering requires, therefore, that 
these phases (pj are adequately restricted, i.e., to an interval of the order around 
a common mean value. A random phase corresponds to a phase variance A0ra = 
tt/Vs. (We use the phase deflnition by Barnett et al. [Q; the numerical factors 
are, however, without influence on the critical part of the entropy.) Restricting 
the phase (pj of the jth atom to the interval (^ corresponds to a reduction factor 
A0k/A0ra- For all atoms this leads to the entropy contribution 

5p.o. = fcB In n (^) ' ^ 5^ nk In (n^) {n^ » 1). (8) 
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This expression appears to be the most simple and plausible way of accounting for 
the phase restrictions. 

Phase ordering requires that the considered phases are reasonably well defined, 
i.e., A0k <^ 1 or fT-k ^ 1. The condition ^ 1 has been used in the last step in 
Eq. (H). We note that the critical terms in the expectation value of the entropy 
are solely due to the contributions from fik 1. The expectation value of Sp ^. 
will exhibit a logarithmic singularity (Sec. |3.3|) , i.e., the phase ordering entropy is 
the decisive contribution in our model. 

For nk 3> 1, expression (|^) yields S'ibg ~ ^felnnk. From this and Eq. (||) 
it follows that the low k contributions to the sum S'ibg + Sp,o. are negative which 
might seem disturbing. The entropy must, however, not be attributed separately 
to each k level; it rather results from the distribution of all particles onto the 
available levels. Consider, for example, an isolated level with the single particle 
energy ei that is occupied by rii atoms. Then we may attribute the energy eiUi 
to this level but we cannot say that In rii is the entropy of this (isolated) level. 
Alternatively one may argue as follows: The low k contributions to S'ibg (or >S'i'bg) 
are small and noncritical. Consequently, any entropy expression that leads to the 
experimentally observed critical behavior (i.e., to Sexp oc — tln|t|, where t is the 
relative temperature) must override the low k contribution of S'ibg (similarly as it 
is done by S'p.o.). 

3 Thermodynamic entropy 
3.1 Statistical assumptions 

In order to determine the thermodynamic entropy we need the expectation values 
(rik) of the occupation numbers. For the wave function (^ with the single particle 
functions (|lD and with hard sphere Jastrow factors, Chester obtained the IBG 
values, i.e., (nk) = (?T-k)iBG- This means that the IBG occupation numbers remain 
valid in spite of hard core interactions. The considered phase ordering correlations 
are small in the sense that \Ecorr\/N <^ k^Tx- Therefore, it appears to be a sensible 
approximation to use the IBG occupation numbers for the phase ordered states, 
too. The phenomenological assumption 

(^k) = (?^k)lBG (9) 

is the basis of our description of the phase transition. It implies that we consider an 
almost ideal Bose gas model. This classification refers to the statistical assumption 
(y) but not to an assumption about weak interactions. 



7 



The thermodynamic entropy is obtained by 

^(T,\/,iV) = (5K)) = 5(K)), (10) 

where 

^(nic) = ^{bg + -^p.o.. (11) 

The contributions S{-qq and Sp,o. are defined in Eqs. (H) and (H), respectively. 
For the relevant states (rik ^ 1) the quantum mechanical variances Arik/'^k ~ 
1 / ^ 1 are small compared to the statistical variances. Therefore, we omit a 
distinction between (fik) and (rik). 

The average occupation numbers of the IBG are 

1 1 
^'"'^ " exp [{e^- fi)/kBT]-l " exp {x^ + r^) - 1 ' ^^^^ 

Here = f^k"^ /2m is the single particle energy and /i is the chemical potential. 
We use the dimensionless quantities x and r, 

2 A^fc^ 12/^ /TON 

a; = — = — ; — and r = — - — — ■ , (16) 
kuT Att knT ' ^ ^ 



where A = 2TTh/ ^/2'Kmk-BT denotes the thermal wave length. The transition tem- 
perature Tx is given by the condition 

\{Tx) = [va^/2)f\ (14) 

where C(3/2) ~ 2.6124 denotes Riemann's zeta function and v = V/N is the 
volume per particle. In the following we use the relative temperature 

t = tiT,v) = ^^^. (15) 

The chemical potential vanishes at the transition point. For |t| <^ 1 the particle 
number condition = XK'^k) yields 

^^'^ = 'lk^ = [a'\t\+b't^ + ... (t<0). ^''^ 



The coefficients are known for the IBG, in particular a = (3/4)(^(3/2)/y^ and 
a' = b' = . . . = 0. As a generalization of the IBG we admit r = a' \t\ + . . . with 
a' > for t < 0. This generalization does not change the character of the transition 
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(neither the point of transition nor the critical exponents). Using the second hne 
in Eq. (0) and the particle number condition we obtain the condensate fraction 

N V2 C(3/2)y' ' ■ 

Here no(t) is the expectation value of the occupation number of the lowest energy 
level. 

There are a number of reasons for generalizing the IBG by admitting t = a' \t\ + 
. . . with a coefficient a' > in the expansion (pIBI): Formally, this generalization 
leads to a more symmetric form of this expansion. Physically, a positive value of a' 
implies a temperature dependent energy gap A ^ kBTa'H"^ between the condensate 
and the noncondensed particles and implies that the lowest level is more rapidly 
occupied than in the IBG. Such a behavior leads to a better agreement between 
the calculated and the experimental temperature dependence of various quantities, 
in particular of the superfluid density [^. As a last argument we mention that 
a' > removes the unphysical divergence {l/k"^ for /c — 0) of the static 
structure function of the IBG for t < 0. Our main result will not depend on the 
parameter a'. 

The well-known expectation value of the IBG entropy (^ is given by 

Sibg{T,V,N) _5v _ , _2 



NkB 2 A3 

where 



95/2{t) + T , 



oo / 2 

exp(— nr 



n=l n=\ 



defines Riemann's generalized zeta function. Usually z = exp(/5/i) = exp(— r^) is 
taken as the argument of this function. We prefer the argument r because of its 
close relation (r ~ \t\) to the relative temperature. 

For the IBG entropy one has to use r = in Eq. ( [T8| ) for t < 0. The modification 
T = a'\t\ with a' > leads to a specific heat that is continuous at Tx (as in the 
pure IBG) but that falls off more rapidly for t < 0. 



3.2 Mean phase contribution 

Because of the common mean phases the IBG expression is replaced by S{^q, 
Eq. (H). The expectation value of the entropy S[^q is 

S{^^iT,V,N) _ 1^ 



k 



2(nk)) ln(l + 2(nk)) - 2{n^) ln(2(nk)) • (20) 



9 



For a first, crude estimate we take into account only those contributions that come 
from (nk) > 1. This yields S'ibg ~ ^BZlkl^W) and S{^q ~ (^b/2) Xlk ^^(^k), 
i.e., ^{bg(T, V, N) ~ 5ibg(T, V, N)/2. 

The actual results for S'ibg and S'j'^q near the transition point are 

Sibg{T,V,N) ^ 1.28iVA;B + + , (21) 

S[^GiT,V,N) ^ 0MNkB(^l + ^t + O{t^)y (22) 

Eq. (pTD follows from the expression ([l^). The value of S[-qq at the lambda point 
has been determined by a numerical integration of the r.h.s. of Eq. (pO|). The 
next term in the expansion is due to the prefactor f/A^ oc T'^^^ in Eq. (|1^). This 
prefactor originates from the replacement of the momentum sum by an integral 
over dimensionless variables; it occurs in both entropy expressions in the same 
way. The higher terms O(t^) in Eqs. (^) and (^) do not agree; moreover, the 
coefficients of these t"^ terms are different for t > and t < 0. 



3.3 Phase variance contribution 

The expectation value of the phase ordering entropy S'p.o., Eq. (|), reads 

= Hn^) (1^1 < O-l)- (23) 

^ k 

For the phase ordering we required A0k ^ 1 in Eq. (||), lets say A0k < 0.1 or 
(rik) > 10^. Because of (nk) ~ l/(x^ + r^) this implies x < 0.1 and r < 0.1 or 
\t\ < 0.1. The condition \t\ < 0.1 indicates the temperature range in which phase 
ordering and the corresponding entropy contribution are expected to be relevant. 
For evaluating the entropy (^) the momentum sum is replaced by an integral. 

This step has to be accompanied by a separate consideration of a potential con- 
densate fraction; this will be done below in Eq. (|28|) . 

For a subvolume of size Vq the r.h.s. of Eq. (plf ) would be Vo/(27r)^ Jo^dk Ank'^. 
The subsequent summation over all subvolumes yields the factor V/Vq] this sum- 
mation has been taken into account in the prefactor on the r.h.s. of Eq. (123). In 
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Appendix ^ we discuss the validity of the step (0) with respect to finite spacing 

1 /3 

Ak = tc/Vq of the momentum values. 

We insert the replacement (|2^) , the dimensionless quantities (0) and the oc- 
cupation numbers (|T2|) into Eq. (^): 

gp.o.- 2 " r,,,-M^-^}f^::)-^i^^^j(rf. (25) 



NkB V^X^Jo exp(x2 + r2)-l ^FA^ 

The critical r dependence of the integral J(r) can be determined analytically (App. 
0), 

J(r) = J(0) - Trrlnr + 7rr[l - ln(2)] + 0{t^). (26) 
Using this, r = (3 C(3/2)/4v/?) \t\ + 0{t'^), and f/A^ = l/C(3/2) + 0(t) we obtain 

^p.o,(T,y,iV) = -^tln|t| ±... (|t|<0.1) (27) 

for the leading singularity of the phase ordering entropy for t > 0. As we will see, 
this result holds for t < 0, too. 

For t < the prescription ( p^ ) has to be accompanied by the replacement 

K) (nk)+no(t)(5(k). (28) 

A phase variance A0o ~ ~ stands for a macroscopic phase coherence. 

This is appropriate for a superfiuid fiow but not for the considered local phase 
ordering. For the phase variances (0) we use, therefore, the continuous part of 
the occupation numbers (rik) only. The contribution 3^°^'^ in Eq. (^) due to the 
condensed particles reads then 

This can be evaluated by using Eq. (|T^) for nQ(t)/N and Eq. (|T^) for r(t). 

For t < the contributions ( pSf ) and ( p9D have to be added. The a' term of 
no{t)/N = [3/2 + 2v^a7C(3/3)]|t| + . . . in Eq. (|9D cancels exactly the logarithmic 
term in Eq. (pSj). The surviving leading term comes from nQ{t)/N = {3/2)\t\ + . . . 
and ln[exp(r^) — 1] = 2 In |t| + . . . yielding again Eq. (^). 

3.4 Specific heat 

The specific heat per particle may be written as 

= 1^^^ ^ r -Aln|t|+i? + ... (t>0) 

N\dtjy^ \ -A'ln\t\+B' + ... it<0). ^^''^ 
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From Eq. (^) we obtain the theoretical amphtudes 



^theor — ^theor " ~ ^B- 



3 
2 



(31) 



If the experimental data for cy are fitted by the form one finds A^^p ^ 0.63 
and A'^^p ^ 0.59 fce; these values are taken from Eq. (4.2) (for low pressure) of 
Ref. ||2^. It is remarkable that we obtain a parameter free result and a sensible 
size for these amplitudes. The experimental ratio Aexp/^exp deviates from 1 by a 
few percent. In our model, the contribution (^) of the noncondensed particles to 
Sp,o. has different signs above and below the lambda point. A correction in this 
contribution could, therefore, lead to a deviation from AtheorMlheor — ^■ 

Experimentally the specific heat cp at constant pressure is more readily accessi- 



ble. It is this quantity that has been measured ||2^, |2^, ^ |2^ over several decades 
of the relative temperature, eventually down to |t| = 2 x 10"^ in a microgravity 
experiment [27|. In the context of these measurements and their analyses we note 
the following points: 

1. Theoretically, a logarithmic behavior of one of the two quantities, cp or cy, 
implies p8|, EOl an analogous behavior of the other quantity in the experimen- 



tally accessible range. The statement |Q that only one of these quantities 
may diverge for |t| — applies to very low (experimentally inaccessible) |t| 
values only |^ . 

The experimental values for cp may be fitted by the form (|30D, too. The 
corresponding coefficients (as given by Eq. (47) of Ref. [^) differ from that 
for Cy by about 5 to 10%. For comparing our result (|30| ) with the experiment 
we may, therefore, also consider the experimental cp. 



2. If the critical exponents a and a' are used as fit parameters one finds 



26, 27 1 small negative values. The following points indicate that this does 
not rule out a logarithmic singularity (corresponding to a — and a' 0) 
for the leading term: 



Arp |]HT| fitted the data of many experiments in order to obtain an equa- 
tion of state for helium. He considered specifically the question of the 
critical exponent a (Sec. 13.1.1 of Ref. ||3l|) and found no statistical sig- 
nificance for a deviation from a leading logarithmic form. Consequently, 
he adopted the value a = a' = (standing for a logarithmic function) 
in his expression for the specific heat cy- 
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(b) By the scaling law a' = 2 — Su the critical exponent for the specific 
heat for t < is related to that of the superfiuid density ~ I^T- The 
deviations from a' = and z/ = 2/3 found in standard fits are consistent 



| 27| with this scaling law. 



In Ref. ||2T[] we presented an alternative fit formula for the superfiuid 



density in which the leading exponent u equals exactly 2/3. This fit 
formula reproduced the data better than standard fit formulas. This 
means that the value z/ = 2/3 (corresponding to a' = 0) is not ruled out 
by the experiment, and that the value of the leading exponent (found 
in a fit) might depend on the higher order terms (used in the fit). 

In view of these points, we consider it to be an open question whether the true 
behavior deviates from a logarithmic one in a way that cannot be accounted for 
by higher order terms. 

We turn now to the coefficients B and B' of the specific heat (|30D. From S{^q, 
Eq. (p2D , we get a contribution of 1.44^3 for both, -Btheor and -B^hcor- -^l- (13) 
we add the terms that are linear in t; these terms follow from Eqs. ( p5D and (|29|). 
Including all contributions we obtain 

btheor ~ -0.52 kB, 5^heor = "S-T fee- (32) 

For -B{heor used the parameter value a' ^ 3 that is (rather uniquely) deter- 
mined by adjusting the model expression for the superfiuid density to the 
experimental temperature dependence. 

The value for -Btheor compares reasonably well with the experimental value 
-Bexp ~ — 0.84/cb (from Eq. (4.5) of Ref. ||2^ for low pressure). 



There is a large discrepancy between -Btheor experimental value p3 

-Bgxp ~ 2.0 ks- In the IBG, the coefficient of the t"^ term in the free energy Fibg 
is the same above and below the transition. This means that the IBG does not 
reproduce the jump of the specific heat that normally accompanies the occurrence 
of a macroscopic order parameter field. For a crude estimate of this missing jump 
we consider the Landau free energy FLandau/ (NkBTx) = rt\ilj\'^+u\ip\'^. This Landau 
model yields a jump of the specific heat, Acy/fcB = r^/2u. The identification 
(|^/,|2) = no{t)/N relates the coefficient in = {r/2u)\t\ to that in Eq. (0). 

This connection leads to Acy = r [3/2 + 2y/7ia' /({3/2)]. Using again a' ^ 3 yields 
Acv ~ 5.6 r. For the parameter r in the expression for -FLandau/(-?\^^BFA) we may 
expect r ~ 0{1). The resulting Acy ~ 5.6 k^ could close the gap between B'^^^^^ 
of Eq. dH) andS^p. 
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3.5 Entropy at the lambda point 



As a last point we consider the value of the entropy at the lambda point. Evaluating 
the phase ordering entropy (p5[) as it stands yields 



5p.o.(Ta) = ■j'^^y. . iVfcB ~ -0.51 iV^B. (33) 

This result contains, however, contributions from momenta for which the condition 
(rik) ^ 1 is not satisfied. 

We estimate the value Sp^ {Tx) that comes from the contributions with (nk) ^ 
1 alone. For small and positive t values we may assume the form Sp^ojN ^ 
S^l {Tx)/N - Aheortln(t) + (5p.o. + Aheor)t. FoUowiug Eq. (^) we argued that 
the phase ordering entropy should be relevant in the range \t\ < 0.1 only. This 
implies 5p.o.(t = 0.1) ^ leading to S^^'HTx) ^ -0.3 A^fce. 

We conclude that the result (|33D might be too large by roughly a factor of 
2. One may, therefore, expect that the contribution Sp,o,{Tx) accounts for the 
difference between S[^q ~ 0.96A^A;b and the experimental value |^ ScxpiTx) ~ 
0.76 NkB. 



4 Concluding remarks 

We summarize the main features of the presented phenomenological many-body 
model: 

1. Following Chester p we use the IBG many-body states multiphed by Jastrow 
factors. 

2. We assume that local phase correlations are the relevant correlations near 
the lambda point. These correlations are specified by the single particle 
functions (0). 

3. The local phase ordering can be described by coherent states. The phase 
variances of these coherent states are the basis for the statistical counting of 
the phase restrictions. This leads to the expression (H) for the phase ordering 
entropy. 

4. For evaluating the phase ordering entropy we use the IBG expression for the 
average occupation numbers. 
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It is clear that global phase ordering plays a decisive role in liquid helium below the 
transition point. This makes our assumption that local phase ordering are relevant 
correlations near the lambda point to some extent plausible. We have proposed 
a specific kind of the local phase ordering in the framework of an almost ideal 
Bose gas model. The use of the IBG occupation numbers is a phenomenological 
assumption. 

Our model yields a logarithmic singularity of the specific heat in a straightfor- 
ward and rather simple way: The expression (P) for the phase ordering entropy 
follows from the phase variances (0). Evaluating this expression with the IBG 
occupation number yields the logarithmic singularity ( p7[ ) (immediately for t > 0, 
and after a slight generalization of the occupation numbers for t < 0). 

We summarize the main results of our almost ideal Bose gas model: 

1. The model yields a realistic expression for the specific heat of liquid he- 
lium at the A transition. For the strength of the logarithmic singularity the 
remarkable result A = A' = 3 is obtained. 

2. The model retains the essential features of the IBG. Therefore, it strengthens 
the suggested relation [|I], 0] between the Bose-Einstein condensation and the 
A transition. 

3. The model offers an intriguing picture for the relevant correlations near the 
A point. This picture provides also perspectives beyond the specific heat 
(some of which are indicated in App. 0). 

Various other quantities may be calculated in the framework of our almost ideal 
Bose gas model: In an earlier paper 0] we calculated the energy as a function of 
the temperature; the present derivation of the specific heat appears to be more 
direct. The critical exponent (3 = 1/2 of the condensate fraction may be reconciled 
with the actual behavior of the superfluid density by assuming that noncondensed 
particles move coherently with the condensate |2l[]. This idea leads to observable 
consequences that have been discussed in Ref. [BB . 
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A Coherent states 



In order to justify the use of relation we establish a connection between the 
phase ordered single particle functions @ and coherent states. Our procedure is 
analogous to Anderson's construction of coherent states for the condensate. 
The many-body wave function (^) can be written in the form 



no no+ni 



^ = SFl[M^j) n (34) 

j=l j=rao+l 

The single particle functions (phase ordered or not) are denoted by ipi, where 
i = 0, 1,2, . . . follows the energy sequence. Without the Jastrow factor the off- 
diagonal single particle density reads 

^(r,r') = <^^|^/'+(r')V^(r)|^^ =r2o<Po(r)Vo(r')+^iV^i(r)Vi(r) + ... , (35) 

where and ip are the particle creation and annihilation operator, respectively. 
Including the Jastrow factor the condensate contribution in Eq. ( ^5]) is depleted 
1 3^ , p5| from uq to a lower value Uc- Similar effects are to be expected for next 
low-momentum contributions. 

All terms in Eq. (BH) besides the first term vanish for |r — r'| oo. This 



off- diagonal long range order 0] is considered |^ as the decisive criterium for 
superfluidity. 

Anderson related the off-diagonal long range order to the mean fleld aspect 
by constructing localized coherent states. The coherent states |coh) are deflned 
such that (coh|-?/;|coh) oc exp(i0), where is the phase of the condensate wave 
function ipQ oc exp(i0). Accordingly, we will construct coherent states |coh) for 
which (coh|'?/'|coh) is proportional to the phase ordered single particle function (Q). 
Anderson's aim was to show that the phase of the condensate and the correspond- 
ing particle number may be treated as macroscopic variables. Our aim is to justify 
Eq. (^ by relating the wave functions (||) to coherent states. 

Anderson |^ introduced finite volumes AV in which the phase of the con- 
densate wave function yjo is approximately constant. We consider phase ordered 
single particle functions that are restricted to finite subvolumes of size Vq. In the 
following we identify AV with Vq and treat both cases simultaneously. 

The single particle functions (pk shall form an orthonormal set in one subvolume 
of size Vq. Then we may write down the following relations between the field 
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operator ip{r) and the annihilation operator Ck: 

^(r) = X]</'k(r)ck, (36) 

k 

Ck = / d'ripl{r)^{v). (37) 

JVo 



'Vo 

A state with ra^ particles of momentum k is given by 



1^) = l^k) = — . c,|; I vac), (3^ 



where |vac) is the vacuum state. In the following we restrict ourselves to one 



specific momentum and omit the index k. The coherent state |3^ is constructed 
as 

CO ^ 

|coh) = exp(-|z|V2)^^|n). (39) 



n=0 



This state depends on a complex number z. Because of n = (coh|cj^Ck|coh) = 
we may set 

z = ^Jfi exp(i0). (40) 
For the state (H) with Eqs. (||) and (ID we find 



(coh|'i/'(r)|coh) = \pn exp(i0) V9k(r) (41) 

for the expectation value of the field operator. 
Using 

- 1 (r ^ ^o) , ^^^^ 

for the lowest state, Eq. (pj) becomes (coh|'?/'|coh) = \/nlVQ exp(i0). So far we 
have reproduced Anderson's construction. 

We are now going to connect the phase ordered single particle functions (|^) 
with coherent states. For this purpose we introduce the single particle functions 

V9k,+ (r) oc exp [i (A;2l/ + fcsz)] exp(+igx) (43) 
V9k_(r) oc exp [i (A;2l/ + A;3-2)] exp(-igx) (44) 

These functions shall be orthonormalized in the considered subvolume. Analo- 
gously to Eqs. (|38| ) with (^) we introduce the n-particle states |n)+ and |n)_ that 
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correspond to (y9k,+ and v^k,-, respectively. The coherent state 



Icoh) = !^^pH^ z^\nU-_z_-\n). ^^^^ 



n=0 

with z = v^exp(i(/)) yields then 

(coh|'?/'(r)|coh) = v^V'k" (^)- (46) 
The phase variance of the coherent state p5|) is given by 



A0k = ;^7=. (47) 

Here we have attached the index k again; the above discussion referred to a single k 
value. Eq. (^) is a well-known result ||2^ for a coherent state of the standard form 
(p^). For the state ( ^5]) we may write |coh) = (|coh)+ — |coh)_)/-\/2 in an obvious 
notation. The phase operator acts separately within each of the two sets of states, 
{|?T,)+} or {|n)_}. Moreover, the states of one of these sets are orthogonal to that 
of the other set (because of the functions (|43|) and (|^) are mutually orthogonal). 
Therefore, the phase variance of the state ( ^5l) is one half of the sum of the phase 
variances of the states |coh)_|_ and |coh)_ that are of the standard form (p9|). This 
means that we obtain the standard result (^) for our somewhat special coherent 
state (P3), too. 



The reason for establishing the relation between the phase ordered single par- 
ticle functions (^) and coherent states ( ^3[ ) is the necessity to justify the relation 
(^) or, equivalently, Eq. (0). The relation is the basis for the phase ordering 
entropy (^). 

Instead of |coh) = (|coh)_|_ — |coh)_)/A/2 we may also consider the coherent 
state (|coh)_,_ -|- |coh)_)/\/2. This orthogonal state corresponds to a cosinus in Eq. 
(0) instead of a sinus. The considered correlation effect requires that only one of 
these sets of states is occupied. This 2 : 1 reduction has been discussed in Section 



2.2 



With respect to the condition rik ^ 1 we note that the average occupation 
numbers (rik), Eq. (|12[), do not depend on the size Vq of the subvolumes. The 
finite size implies, however, a finite spacing of the momentum values (App. P) 
that takes care of the finite number Nq = Vq/v of atoms in one subvolume. 
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B Phase coherence volume 



We start by explaining why finite phase correlations require that the single particle 
functions (@) have to be localized within finite volumes. Then we show that the 
lowest single particle functions may extend over several of these subvolumes. This 
leads to the notion of a phase coherence volume for which we obtain V^oh ~ ^o/|^P, 
where Vq is the size of one subvolume. We discuss this result in a number of points. 

Two phase ordered states with momentum k and k' are correlated if ki = k[. 
The restriction ki = k[ cancels one momentum sum in the expression for the 
energy. A momentum sum fdk/Ak is proportional to = V^^^/27i 

or to iV^/^. The correlation energy -Ecorr is, therefore, proportional to N"^/^ only. 
This is consistent with the observation that a phase ordering of the kind (|]) can 
be obtained by physical boundary conditions at the wall of the considered volume. 
For an infinite volume such a surface effect vanishes (Ecorr/N oc — > 0). A 
finite correlation effect may, however, be obtained by dividing the total volume V 
into V/Vq subvolumes of size Vq and by requiring physical boundary conditions at 
the walls of these subvolumes. For the whole system (including a factor V/Vq for 
summing over all subvolumes) we obtain then a result of the form 

^corr Wo ^^^^ 



N 



where < is some average strength of the attractive part of the interaction (for 
a more detailed discussion see Ref. [^]). This result shows that a finite correlation 
effect requires a finite size Vq = iVg v of the subvolumes. 

Finite volumes Vq imply a finite spacing Ak = tt/Vq^^ of the momentum values. 
Alternatively one may start from a finite spacing and admit only the momentum 
values 

qn = qo + n-Ak, (n = 0, 1, 2, . . .), (49) 

in the single particle functions (^. The spacing ( ^9]) leads to a finite correlation 
energy even if the single particle functions are not localized. In a macroscopic sys- 
tem the possible q values are, however, dense, and the entropy drives the particles 
into the occupation of all available states. This is the reason why the finite spacing 
(^) can be realized only for single particle functions localized within subvolumes 
of the size Vq = (vr/AA;)^. In this sense, a finite spacing implies finite volumes. 

This discussion shows that we may either start from finite volumes or from a 
finite Ak value. There is, however, the following difference: Starting from the finite 
spacing (^9]) we may admit a go value below Ak without damaging the correlation 
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energy. A single particle function with go < can, however, not be localized 
within Vq. 

In view of this observation we consider the following modified picture. Only the 
single particle functions (0) with g„>i are localized within Vq, the single particle 
function (fQ with go may have a larger extension. We present an estimate for the 
volume V^oh of the lowest single particle function ipQ-. Let V^oh be some multiple of 
Vo, that means Koh = WVq. In this case, there are W single particle states with 
g < Ak in the volume V^oh out of which only one (the one with go) is occupied. 
A redistribution of tiq atoms over these W states would increase the entropy by 

As ~ /cb ln(r2o)'^- At the same time these atoms would loose their correlation 

1 /s 

energy, Accorr ~ —noWo/NQ . The stability condition Ta As < Accorr yields the 
upper bound W < [—WQ/{NQ^^kBTx)]no/\nnQ ~ tiq. Using Hq = (rik^o) = ~ 
leads to 

Koh(t) = W^^~^. (50) 
We discuss this result in a number of points: 

1. The lowest single particle functions v?o extend over volumes of the size Vcoh- 
Within a volume V^oh the function ipo defines the direction of the phase 
ordering (x direction in Eq. (0)) and the phase 0. Therefore, Vcoh{t) of Eq. 
( pOD constitutes a phase coherence volume. 

2. Strictly finite volumes for all particles would mean that the lower bound of 
the integral J{t) in Eq. (^) would be Ax instead of zero. This would imply 
a cut of the logarithmic singularity. The replacement (|2^) may, however, be 
justified by a lower bound go ~ 1 /Koh and an average over somewhat different 
go and Ak values in different parts of the macroscopic system. A closer 
examination of this point might lead to a modification of the logarithmic 
singularity at very low \t\ values. 

3. In the estimate leading to Eq. (0) we used the continuous part of the occupa- 
tion number (nk^o) also for t < 0. Using instead the condensate occupation 
number no(t) leads to an infinite volume, V^^^ = oo. This is the adequate 
phase coherence volume for the potentially macroscopic range of a superfluid 
flow. 

4. Different phase coherence volumes of size Vcoh within the macroscopic sys- 
tem will, in general, correspond to different phase directions (x direction in 
Eq. (j^)). A specific phase direction implies a small anisotropy of the static 
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structure function. In principle, the volumes V"coh and their increasing size 



for i ^ should, therefore, be observable. 



5. When approaching the lambda point from above the phase coherence vol- 
umes grow according to Eq. (|SD|). At the lambda point the coherence volume 
becomes infinite. This picture has some similarity with a ferromagnetic sys- 
tem. 



Below the lambda point the finite volume ( pOD refers to the phase ordering 
of real single particle functions of the form (|^). At the same time there 
is a potentially infinite range of phase coherence of a superfiuid current. 
The corresponding phase coherence volume V^^-^ — ^ refers to the complex 
condensate wave v^o oc exp[i0(r)]. 

6. Adjusting the correlation energy (ffSp to the experimental strength of the 
logarithmic singularity yields |Q 50. The phase (p should be approxi- 

mately constant within Vq (see App. 0). It appears, therefore, tempting (but 
also speculative) to set |V0|max ~ Vq and to use the relation Ug = hV(j)/m 
for the velocity of a potential superfiuid flow. The maximum (or critical) ve- 
locity for such a flow would then be 

^|V(/>|max h hi m 

Vcr\t = YE = 175^7^ ~ 0-8 — • (51) 

m mV^ mv^i'^ ]\f^''^ s 

In this way the finite size Vq renders a possible connection between the "nat- 
ural" velocity scale h/{mv^^^) and realistic values for the critical velocity. 
Using Eq. (|50D we obtain t>crit ~ |t|^^/^m/s for the critical velocity near the 
transition point. 

C Evaluation of J(r) 

We determine the critical (r ^ 1) behavior of the integral 

/•°° , 2 ln[exp(x2 + r^) - 11 
Jr=/ dxx' ^ ^\ ' r 52 
Jo exp(x^ + r^) - 1 

Writing ln[exp(x^ + r^) — 1] = ln[l — exp(— — r^)] + + we obtain 



~ x^ 



ln[l - exp(-x2 - r^)] p x^ (x^ + r^) 



J(r)= / dx ( 2 ^ 2\ ^ + / / 2 I ^■ (53) 

exp(x^ + r"=) — 1 Jq exp(x"^ + r^j — 1 
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The second integral yields a result that is of the form 

i?o = const. + C(r2). (54) 

In the following Rq stands for any expression of this kind. By O(t^) we mean 
terms that are proportional to or to higher powers of r. 

In the first integral in Eq. (^) we use the following expansions into powers of 
y = x'^ + or exp(— y): 

f ln(y) Hy) + 1 , y (2 Hy) + 7) ^ 

ln[l - exp(-y)] ^\ y 2 + 24 ' ' 

exp(y)-l I 3exp(-3i/) llexp(-4y) 

1^ - exp(-2y) ± . . . {y>l). 

(55) 

We divide the integration into one part from to 1 and another part from 1 to 
oo, and insert the appropriate expansion (0). All terms except the one with 
\n{y)/y yield contributions of the form (Q). We evaluate the remaining integral 
over ln(?/) /y = ln(y) — ln(?/) /y, 

J(r) = f'dx Hx' + r2) - Cdx ^""^f ^ + 
Jo Jo + 

, r , ln(x2 + r2) 

The first integral yielded vrr + /2o (see number 2733 in Ref. [^). Because of 
Jj^ (ix ln(x^ + r^)/(x^ + r^) = Rq the upper limit of the remaining integral could 
be set equal to infinity. In this integral we substitute x = tz: 

J(r) = nr — T dz — r / h ito 



l + z' Jo 1 + z' 

= J(0) -7rrlnr + 7rr[l-ln(2)] + 0(r2). (57) 

The first integral is elementary, the second one may be found under number 4.295 
in Ref. |38]. A numerical integration yields J(0) ~ —1.183. 
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